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Abstract 



Quantum Heisenberg ferromagnets with long-range interactions in low di- 
mensions are investigated by means of the Green's function method. The 
model Hamiltonian is given by H = —(Jo/2) J2i^j r ij P ^i ' &j ~ It 
is shown that there exists a finite-temperature phase transition in the region 
d < p < 2d for the (i-dimensional case and that no transitions at any finite 
temperature exist for p > 2d; the critical temperature is also estimated. We 
study the magnetic properties of this model. We calculate the critical expo- 
nents' dependence on p; these exponents also satisfy a scaling relation. Some 
of the results were also found using the modified spin-wave theory and are in 

remarkable agreement with each other. 
PACS numbers : 75.10.Jm, 05.70.-a, 75.40.Cx 
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I. INTRODUCTION 



Models with long-range interactions (LRI's) have attracted much attention since the 
exact solution of the 5=1/2 Heisenberg model with interactions proportional to the inverse 
squared distance between sites (HS model) was found independently by Haldane |T[ and 
Shastry 0. In the case of models with interactions decaying as r~ 2 we have made use of the 
mathematically simple structure of those models to obtain many exact results. We cannot, 
however, apply the same method used for the r~ 2 case to the case when the interactions 
have the general form r~ p , because the mathematical structure of this general case is more 
complicated. Among ferromagnetic models with r~ p interactions, the Ising model, classical 
Heisenberg model, n-vector model, and spherical model have already investigated. Regarding 
the Ising chain |3|-|5| , the existence of a finite-temperature phase transition (FTPT) when 1 < 
p < 2 and its absence when p > 2 are known. The region where there exists a FTPT in the 
one- and two-dimensional case of the classical Heisenberg model and the spherical model 

is known to be d < p < 2d, where d denotes the lattice dimensionality which is taken to be 

1 or 2. The thermodynamic properties of the <i-dimensional spherical model were discussed 
0. Critical exponents of the n-vector model were obtained by means of the renormalization 
group approach A Monte Carlo simulation of the classical Heisenberg model was 
recently reported [|T^]. However, less is known about quantum Heisenberg ferromagnets 



with LRI's. The absence of a FTPT in the region p > 2d [T3J and upper bounds for the 
correlation functions [12 1 are known. In this quantum case, however, little is known about 
the thermodynamic properties and critical phenomena. 

Recently, we have reported the low-temperature properties of the one- and two- 
dimensional quantum Heisenberg ferromagnets with LRI's decaying as r~ p , using the mod- 
ified spin- wave (MSW) theory [jT3,p~4 1 and the Schwinger-boson mean-field theory |15 |. In 
those papers, the critical temperature of the FTPT, the temperature-dependence of the 
susceptibility and specific heat were obtained. In particular, from both of these two approx- 
imate methods, results in the one-dimensional case for p = 2 are in good agreement with the 
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exact solution of the HS model. We cannot, however, obtain the properties of the ordered 
phase of this model. The reason is that these methods are approximations approaching from 
the disordered phase to the critical point. 

In this paper we will apply the Green's function method (GFM) with Tyablikov's de- 
coupling [IB| to the above model in an external magnetic field. Its Hamiltonian is given 



by 



w = 4 E E J(p)Ss ■ s j+p -hy,s*. (i.i) 



2 j P 3 

Here, the strength of the coupling is defined by 

J» <' = <» (1 . 2) 
->0 ° |p| p (otherwise), 

where N is the number of sites and Jo means the nearest-neighbor interaction, which is 
positive. The condition p > d is necessary for sensible thermodynamics. Spin operators 
are taken to have the usual properties: Sj ■ Sj = S(S + 1), [Sj,S%] = iSj5jk, and so 
on. This approximation was first used by Tyablikov in investigation of quantum spin 
systems with the usual short-range interaction and has subsequently been developed by 
many other researchers ||17|| . In those days, however, the three-dimensional case was of 
primary interest and it was only a few years ago that this approximation was applied to 
low-dimensional systems ||18|| . In three dimensions, many good results can be obtained 
with this approximation. Particularly, the critical temperature estimated from this method 
agrees well with other estimations. In the low- dimensional case, on the other hand, the 
critical temperature is found from the GFM to vanish. Although this result T c = in 



low dimensions does not break the Mermin- Wagner theorem [IP , it was considered that we 
could not obtain reliable information in the low- dimensional cases using this method which 
is valid at temperatures under T c . Considering in this paper that the low-dimensional result 
of T c = and the finite value of T c in three dimensions are good, we will estimate the 
critical temperature of the low-dimensional model ( |1.1|) with FTPT. Using the GFM as an 
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approximation approaching from the ordered phase to the critical temperature, in addition, 
we will also discuss magnetic properties of this model. 

This paper is organized as follows. First, a formulation of the Green's function method 
is explained. Then, the equations which the magnetization should satisfy are derived. Using 
these equations, the critical temperature is estimated. The critical behavior of the magne- 
tization and susceptibility are also obtained. Results are compared with previous results. 

II. METHOD 

Let us consider the retarded double-time Green's function (GF) of operators A and B, 
which is defined by 

«i(0;B(o» = -<^-oa^),s(O]>, (2.i) 

where 9(t - 1') denotes the step function. Here, (6) = Tr[6 exp(-/3'H)]/Tr[exp(-/3 / H)] and 
0(t) = e ,w Oe~ !M , where (3 denotes the inverse of the temperature T. Differentiating this 
GF with respect to time t, we find i^({A(t — t'); B)) to be 

5(t - t')([A(t), B(t')}) + (([A(t - t'),H]; B)), (2.2) 

where we should note that ((A(t); B(t'))) = ((A(t — t'); B)) and 5 denotes the delta function. 
Using an approximation, we can find the GF of interest from its time derivative as shown 
above. We note here that the relation between the spectral intensity of the correlation 
function (B(t')A(t)) and the Fourier transformation of the GF with respect to time is known 
as the spectral theorem ans is given by 

Urn [G AB (E + iA) - G AB (E - iA)} 

= (e? E - l)I(E)/i, (2.3) 

where the Fourier transformation Gab(E) is defined by ((A(t — t'); B)) = 
IZo G AB(E)exp[-iE(t - t')]dE and the spectral intensity I(u) is defined by (B(t')A(t)) 
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Here, in order to apply the above procedure to the Hamiltonian ( |1.1|) , let us consider the 
GF (( J S'j + (t); 5~)). The first term in ( [2.2|) is then given by 2(Sf)5(t)Si m and the second term 
is calculated to be 

H((SHt); S-)) + £ J(p)[((St(t)S^ p {t); S-)) - ((SfflS&M S m ))]. (2.4) 



Now, in order to find the GF from its derivative, we introduce Tyablikov's decoupling |L6 
then we approximate ((S?(t)S+ p (t);S-)) and ((S+(t)Sf +p (t); S~)) by (5f ) ((S+ p (t); S~)) and 
(Si +p )((Si(t); S~)), respectively. We can consider that this approximation is valid while (Sf) 
has a nonzero value. Because the value of magnetization (5f ) is considered to be independent 
of its site I, we can set (Sf) = M for any site I. Thus, we obtain for the equation of motion 
of the GF 

i^«S+(f); S-)) = 2MS(t)8 l m + H((S+(t); S~)) + M£ J(p)[((5+(t); S m )) - ((S+ p (t); S m ))]. 

(2.5) 

We now use the Fourier transformation, defined by ((Si~(t); S~)) = J dco-^ J2k e lk ^~ m ^~' ltJt 
G k (u), in order to solve Eq. Q, giving G k {u) = M/[tt(u — H — e(k))}, e(Jfe) = J M[^ p (0)- 
J p (k)], and j p (k) = I] p e iA >p J(p)/ Jo. From the spectral theorem and the solution of the GF, 
we find the correlation function (S~Sf(t)) to be £ fc e ik(i-m)~i[H+s{k)]t / { e p[H+<k)] _ i}. 
When I = m and t — 0, we have 

s(s + 1) _ ((5 , T) ^ ?C0th (ffiL±iWl). (2 . 6) 

In order to discuss the case of general S, it is necessary to consider the Green's functions 
({S+(t); (S z m rS-)) for n = 0, 1, 2, • • • , 2S - 1. From «S+(t); 5* 5~», we follow the same 
procedure giving 

2wsrsf> , 1 = j. Ecot Jfcf(iV (2 . 7) 



3((5f) 2 ) -M- 5(5 + 1) AT 



k 
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III. RESULTS AND DISCUSSION 



First of all, we estimate the critical temperature (T c ) of spontaneous magnetization (i.e. 
H = 0), which we can obtain from Eqs. ( [2.6| ) and (|2.7|) . We take the limit M — > in Eq. 
to have 



[S(S + l)-((Sff)]J /T o = 2I p (0), (3.1) 

where I p (a) is defined by X^[ a + 7p(0) — 7p(&)]~ 1 - I n the same way, we find that ((5f) 2 ) = 
S(S + l)/3 from Eq. (|2.7|) . Hence, we obtain 



Jo _ 3/ p (0) 



(3.2) 



T c 5(5 + 1) 

In the thermodynamic limit N — ► oo, we can replace the summation (l/N)J2k by the 
integral / 1BZ d d fc/(27r) d , where 1BZ denotes the first Brillouin zone. According to Ref.[T3], we 
can calculate the infrared behavior of the function 7p(fc) and obtain 



7 P (0) - l P {k) oc < 



k 2 (p>d + 2) 

k 2 \n(l/k) (p = d + 2) (3.3) 
k p - d (d<p <d + 2). 

Therefore, the integral in Eq. ( |3.2|) is found to be convergent in the region d < p < 2d and 
divergent for p < 2d. This means that the critical temperature is nonzero when 1 < p < 2 
in one dimension and when 2 < p < 4 in two dimensions. 
For the 5 = 1/2 case, we can calculate T c from 

Jo_f 4(27r)- d d d fc _ r 4(27r)- d d d k 

where v4W=rf d " p 7r d [r(p)]-7sin[7r(p-d)/2], 5«=C(p)> and S| 2 ) = C((p/2) -l)/3((p/2) -1). 
Here, T(p) and ((p) denote the gamma function and Riemann's zeta function, respectively, 
and (3(a) is defined by Y2£=i(— l) n ~ 1 (2n — l)~ a . Numerical results for T c are shown in Fig.|l|. 
These results agree qualitatively with those from the MSW theory [13]. Note especially 



that in two dimensions the discontinuity in T c is obtained at p = 4 and that the critical 



temperature vanishes just at that point. In the one- dimensional case, on the other hand, 
there is no discontinuity in T c . Using Eq. ( |3.2| ), we have for the critical temperature of the 
general-5 1 model in one dimension near p ~ 2 

^ = vr 2 ((5f) 2 )(2 -p) = ^S(S + 1)(2 - p). (3.5) 

We should remember that from the MSW theory the critical temperature in the one- 
dimensional case is found to be n 2 J S 2 (2 — p) near p ~ 2 and that S 2 is the dominant 
term of ((S'f) 2 ) under infinitesimally small H in the large-S" limit. We can, then, consider 
that the result Q3.5p also supports the previous result from MSW theory. 

Next, let us investigate the spontaneous magnetization at temperatures below T c in the 
region d < p < 2d. Henceforth, we concentrate on the S = 1/2 case. We have, then, from 
Eq. (pD 



^EMcoshf 



(3.6) 



2 J 

We use Eq. (|3.6|) in order to find numerically the behavior of the spontaneous magneti- 
zation M. Results are shown in Fig.|2|. As we expected, we can see that the smaller p is, the 
more stable the magnetization is and that the two-dimensional magnetization is more stable 
than the one-dimensional one. At low temperatures {(3 ^> 1), we obtain from Eq. 



1 T(d/p-d) ( 2T 



p—d 



M ~ ^ jjr . (3.7) 

2 d7r(p-d) \J A?) y ' 

Near the critical temperature, on the other hand, the spontaneous magnetization is very 
small. The right-hand side (r.h.s.) of Eq. (|3.6|) is then expanded with respect to M and 
e(= 1 — T/T c ); the dominant terms are found to be | — | + J2k{lp(fy ~ 7p(^)]- m this 
expansion, the second term and the third one determine the critical behavior of the magne- 



tization. The magnetization is then calculated to be i/3T c e/2 JqC(p) f° r the one-dimensional 



case and y 3T c e/4JoC(p/2)/?(p/2) for the two-dimensional case. The magnetization exponent 
is found to be 1/2, which is independent of p. This value is the same as the result from 
a mean-field(MF) approximation. However, the MF theory leads to the FTPT even in the 



case of the low- dimensional quantum Heisenberg ferromagnets with usual nearest-neighbor 



interactions by mistake, in spite of the Mermin- Wagner theorem about that model |I9" |. This 
point is a decisive difference between the GFM and the MF theory. 

Now, we calculate the exponent with respect to the temperature dependence of the 
susceptibility near T c . The susceptibility when the applied field is small is defined by x = 
M/H. The summand in the r.h.s. of Eq. ( |3.6|) is then rewritten as ^-(1 + ^ + • • •), where 
x is defined by (3JqM[-^ + 7 P (0) — j p (k)]/2. Therefore, we obtain for the equation which 
determines the susceptibility at temperatures near and above T c 

< 3 - 8 ' 

In the case when p > 2d, the FTPT does not exist, i.e. T c = 0. For p > 3 in the 
one-dimensional case, the leading-order term on the r.h.s. of Eq. ( |3.8| ) is proportional to 
y/xJo- At low temperatures, hence, we find \ oc (3 2 . The Bethe-ansatz method leads to 
the same divergence of the susceptibility at low temperatures in the case of the model with 



nearest-neighbor interactions (p — > oo) |20|j2T[| ; the present result includes it. In the same 
way, we find in the one- dimensional case 

f(/3J )^ (2<p<3) 
X < (3.9) 
[ exp(/?J 7r 2 /4) (p = 2), 

and we obtain \ x ex P (flJo'K((p/2 — l)(5{p/2 — 1)) for p > 4 in two dimensions. These 
results for \ without the FTPT agree with results from the quadratic theory of the MSW 
approximation. We should note that in the case of the HS model (p=l, d—1, and 5=1/2), 
the susceptibility is exponentially divergent at low temperatures, which is in qualitative 
agreement with its exact solution [ ^2] . 



For d<p<2d where the FTPT exists, Eq. (gj) for S = 1/2 and Eq. (|3J) give for the 



equation which \ & t temperatures near T c should satisfy 

^^=/p(o)-/p(4)- ( 3 - iq ) 

4 X J o 
It should be noted here that the function I p (a) behaves for positive and small a as 
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2d~p 

a p~ d 



[\d<p< 2d) 



I p (0) - I p (a) oc < 



The susceptibility is then found to be 



a ln(l/a) (p = |<f) 
a (d < p < 



(3.11) 



p—d 



(f3 c -f3)—* {U<p<2d) 



i 

I 0c- P 



(3.12) 



(d <p < fef). 



In comparison with the magnetization exponent, which is independent of p and which is 
found to be the MF value, it is remarkable that the result ( |3.12|) for the susceptibility 
exponent is dependent on p. 

Finally, we consider the critical isotherm exponent 5' defined by M oc H 1 ^ 5 ' for small H 



and M just at the critical temperature. From Eq. ( |3.6|) , we obtain for small H and M at T c 

(3.13) 



H, ((3 C J ) 2 



12N 



m 2 EWo)-#)]. 

k 



Therefore, we find from this equation 



M oc < 



2d— p 

Hp 



\d <p < 2d) 



(H\n(l/H)) 1/3 (p=ld) 



(3.14) 



#1/3 



(d < p < \d). 



So far, we have known three kinds of critical exponents in the region where there exists 
the FTPT, the magnetization exponent f3', the susceptibility exponent 7', and the critical 
isotherm exponent 5'. It is remarkable that the scaling relation 7' = f3'{5' — 1) is always 
satisfied in spite that these exponents are dependent on p. The same behaviors of these 
exponents are known in the case of the spherical model [0 and the n-vector model |||9| with 
LRI's, but this is the first report of these exponents in the case of the quantum Heisenberg 
ferromagnets with LRI's. 
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IV. CONCLUSION 



We have studied the quantum Heisenberg ferromagnets with long-range interactions in 
low dimensions within the framework of the Green's function method. It is shown that 
even in low dimensions there exists the finite-temperature phase transition in the region 
d < p < 2d; the critical temperature in the region is estimated. Critical behaviors of this 
model are also calculated and found to be depend on p. The present results strongly support 
the previous results from the modified spin-wave theory, which are qualitatively reliable in 
comparison with the exact results in one dimension when p = 2 and when p — > oo. 
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FIGURES 

FIG. 1. The critical temperature estimated numerically. The one-dimensional case is shown by 
the solid line. The broken line denotes the result in two dimensions. 

FIG. 2. Numerical results for spontaneous magnetization of the 5=1/2 model. One-dimensional 
and two-dimensional results are shown in (a) and (b), respectively. 
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